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Abstract 

The aim of this contribution is to give the explicit formulas for the eigenvectors 
"C^." ' of the transfer-matrix of Baxter-Bazhanov-Stroganov (BBS) model (iV-state spin 

model) with fixed-spin boundary conditions. These formulas are obtained by a 
limiting procedure from the formulas for the eigenvectors of periodic BBS model. 
The latter formulas were derived in the framework of the Sklyanin's method of 
separation of variables. In the case of fixed-spin boundaries the corresponding T — Q 
C/3 . Baxter equations for the functions of separated variables are solved explicitly. As 

£S| \ a particular case we obtain the eigenvectors of the Hamiltonian of Ising-like Zjv 

quantum chain model. 

1 Introduction 

« ■ 

During the last two decades, a considerable progress in application of the Separation of 
Variables (or functional Bethe ansatz) method to a broad class of integrable models of 
quantum chains and statistical physics has been achieved. This progress was initiated by 
the paper [lj of Sklyanin who has proposed a recipe for a Separation of Variables in the 
case of quantum Toda chain where the algebraic Bethe ansatz fails. The idea is to write 
the eigenvectors of periodic problem as a linear combinations of the eigenvectors of an 
auxiliary problem (the auxiliary problem for the periodic Toda chain is the open Toda 
chain) . The next step is to construct recursively the eigenvectors of the m-site auxiliary 
problem in terms of the eigenvectors of the (m — l)-site auxiliary problem (2J [3]. 

In [4], this program was realized for the inhomogeneous periodic Baxter-Bazhanov- 
\ Stroganov (BBS) model (or T^ 2 \t g ) model) El [7], [8] defined in terms of cyclic L-operators 

[HI [7] . It was shown that for every eigenvalue (found from the functional relations) of the 
transfer-matrix there corresponds an explicit formula for the eigenvector. So the problem 
of finding the eigenvectors reduces to the problem of finding the corresponding eigenvalues. 

At special values of parameters, the periodic BBS model can be interpreted as a 
model with fixed-spin boundary conditions [5]. In this paper we derive the eigenvectors 
of this model from the eigenvectors of periodic BBS model by specializing corresponding 
parameters. In the case of fixed-spin boundary conditions, the structure of eigenvalues 
is simple, so the T — Q Baxter equations for the functions of separated variables can be 
solved explicitly in terms of cyclic function w p (j) [10] which is a root of unity analog of 
the g-gamma function. It gives us formulas for the eigenvectors of the BBS model with 
fixed-spin boundary conditions. 
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At the end of the paper we show that the Hamiltonian of the Ising-like Zat quantum 
chain model with fixed boundary spins found by Baxter [6j [11] can be obtained in the 
present framework by further specialization of parameters. It means, in particular, that 
we can find explicit formulas for the eigenvectors and the eigenvalues of the Hamiltonian. 
The formula for the eigenvalues coincides with Baxter's one obtained from the functional 
relation for the transfer-matrix. Explicit formulas for the eigenvectors prove that the 
corresponding eigenvalues enter to the spectrum of the Hamiltonian with the multiplicity 
one. 

2 Lattice and quantum chain formulations of BBS 
model 

Following the notation of a recent paper of Baxter [5], we define the BBS model as a 
statistical model of short-range interacting spins placed at the vertices of a rectangular 
lattice. We label the spin variables s x>y by a pair (x,y) of integers: x = 1, ...,n+ 1, 
and y = 1, . . . ,m. Each spin variable 8 x>y takes N values (N > 2): 0, 1, . . ., N — 1. 
The model shall have Z^-symmetry and we may extend the range of the spins to 
all integers identifying two values if their difference is a multiple of N. The model has a 
chiral restriction on the values of vertically neighboring spins: 

S x ,y - $x,y+l = Or 1 mod iV . (1) 

In the following we will consider the spin variables on two adjacent rows: (k, I) and 
(k, l + l), where I is fixed and k = 1, . . . , n+1. Let us denote s^ = 7^ and s^z+i = 7^. The 
model depends on the parameters t q and a' k , b' k , c' k , d' k , a k , b k , c k , d k , k = 1, 2, . . . , n + 1. 
Each square plaquette of the lattice has the Boltzmann weight (see Fig{T]) 

1 

Wr(73fc-i,7fc;y*-i,y*)= E u m ^-^\-ut q y^~ m ^ x (2) 

m fc _i=0 

where 10 = e 2m / N , and 

*Z(0,0) = 1, ^(0,1) = -c^, iftl,0) = i iftl,l) = ^4, 

°k °k °k 

and expressions for F k {j k — 7^, m k _i) are obtained from F k {p/ k — 7^,, m k ) by substitutions: 

h' r' d' —> n" h" r" d" 

a k ) u k 1 c k 5 u k ^ a k ) u k > c k ) u fc • 

We consider the periodic boundary condition: 7„ +2 = 71, 7^ +2 = 7i, where n + 1 is 
the number of sites on the lattice along the horizontal axis. The transfer-matrix of the 
periodic BBS model is iV n+1 x N n+1 matrix with matrix elements 

n+2 

Wi(7,y) = nWrCTfc-i.TJfeJTfc-i.y*), (3) 

k=2 

labelled by the sets of spin variables 7 = {71, 72, ... , 7n+i} and 7' = {71, 7 2 , • • • , 7^+1} of 
two neighbour rows. 
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Fi gure 1: The triangle with vertices marked by the spin variables 7*— i, 7^_i , mk—1 corresponds to the 
function F^^fc-i - 7^, m fc _r) in {2}; the triangle j k , j' k , m fe _i to F k '(j k -j' k ,m k -i). 

Considering m k , A; = 1, . . . , n + 1, in (j2J) as auxiliary spin variables which take the two 
values and 1, we can rewrite transfer-matrix ([3]) in a vertex formulation associating a 
statistical weight not to the plaquettes but to vertices each of them relating four spins: 
m fc _x, m k , 7 fc , j' k (see Fig. 1). Then the weight associated with the kth vertex is 

4(t g ;m fc -i,m fc ;7fc,7fc) = u m ^- m ^(-ut q y^- m ^x 

x F£(j k - j' k , m k „ x )F' k ( lk - j' k , m k ). (4) 
and the transfer-matrix ([3]) can be rewritten as 

n+2 

W(7,7') = Y[h(t q ;m k _ u m k ;-f k ^' k ). (5) 

mi,..., m n +i fc=2 

For our construction of the BBS model eigenvectors we will use a description of this 
model as a quantum chain model. To the each site k of the quantum chain we associate 
the cyclic L-operator [HI U\ acting in a two-dimensional auxiliary space 

/ 1 + Ax fc v fc , AUfc X (a fe - b k v k ) \ 
L h (X)=[ ), fc = l,2,...,n + l. (6) 

\ u k (c k - d k v k ), Xa k c k + v k b k d k /x k J 

At each site k we define ultra-local Weyl elements u k and obeying the commutation 
rules and normalization 

UjU fc = u fc Uj , VjV k = VkVj , Uj-Vfc = -A' v/,u ; , lo = e 2nl/N , uf = vf = 1. (7) 

In (jSj), A is the spectral parameter and we have five parameters x k , a^, b k , c k , d k per 
site. At each site k we define a iV-dimensional linear space (quantum space) V k with the 
basis \j) k , 7 G Zjv and natural scalar product fc(7'|7)fe = 5y )7 . In 14 the Weyl elements 
u k and Vjt act as: 

^k\l)k = ^\l)k, vjfc|7) fc = |7 + l) fc . (8) 
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The correspondence between the lattice BBS model and its quantum chain analog is 
established through the relation 



h(t q \ m k -i,m k ; 7fc,7fe) = k( r rk\Lk(X) mk _ umk \lk}k ( 9 ) 
and the following connection between the parameters of these models 

d' d" r" fj' n" r' n' d" 

\ , ,+ „ a k a k c k l ,, a k a k c k j a k a k / ln \ 

\ = -U)t q , = 77" "TJJ i a k = — 1 b k =U — —, C k = —, d k = ——. (10) 

u k u k U h U k u k u k u, u. 

A A n A- r\j r\: A A 

We extend the action of the operators u^, to y( n+1 ) = Vi ® V2 <E> ■ ■ ■ <8> V n +i defining 
this action to be trivial in all V s with s ^ fc. The monodromy matrix for the quantum 
chain with n + 1 sites is defined as 

T n+1 (A) = L x (A) L 2 (A) • • • L n+1 (X) = ( ^ n+1 ^\ ^ . (11) 



C*n+i(A) -D„, + i(A) 

The transfer-matrix ((Sj) is obtained taking the trace in the auxiliary space 

Wi(A) = tr T n+1 (X) = A n+1 (X) + D n+ i(X) . (12) 

This quantum chain is integrable because the L-operators (jSJ) are intertwined by the 
twisted 6- vertex i?-matrix 



R(X,u) 



( X-ujv \ 

uj{X-v) A(l-w) 

-w) A-z/ 

V X-uu / 



(13) 



J2(A,i/)4 1) (A)L£°(i/) = L? ) (i/)4 1) (A)i2(A,i/), (14) 

where ^ 1} (A) = L fe (A)®I, ^fV) = Relation flS]) leads to [t n+ i(A), t n+1 (fj,)) = 

0. So t n +i(A) is a generating function for the commuting set of non-local and non- 
hermitian Hamiltonians of the model. It also follows from the intertwining relation ffl~4l) 
that [B n+ i(X), B n+ i(fi)] = and therefore B n+ x(X) is a generating function for another 
commuting set of operators. 

The BBS model with fixed-spin boundary conditions can be obtained using periodic 
inhomogeneous BBS model with n + 1 sites (in each row if we consider the lattice for- 
mulation) if one fixes a' n+1 = d' n+1 = 0. In this case F^ +1 (7 n+ i — j' n+ i,m) = unless 
7n+i = l'n+v It means that t n+1 (A) with such parameters and fixed 7 n+ i can be inter- 
preted as transfer-matrix of n-site BBS model with first and last (nth) spins interacting 
with fixed spin 7 n+ i (fixed-spin boundary conditions). 

Let us look what we have in L-operator formulation. Due to (1101) . the relations a' n+1 = 
d' n+1 = give us x n+1 = b n+1 = d n+1 = and b n+1 d n+1 / x n+1 = 0. So we get 



J n+1 



(A) = ( 1; ) = ( 1 > ) ( 1, Au-> n+1 ) . (15) 

u n+1 c n+ i, Aa n+ iC n+ i 1 \ u ra+1 c n+ i 



Since t n+1 (A) commutes with u n+1 we may restrict the space of states of (n + l)-site 
problem to the states with fixed eigenvalue of u n+ i = u ln+1 in order to obtain ra-site 
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problem. To consider different fixed spins on both boundaries we include one more L- 
operator Lq(A) like (fT5|) but with operator Uo and parameters ao and cq. On the states 
with fixed eigenvalues u n+ i = cj 7n+1 and u = u 10 we have 



L n+1 (A)L (A) = (1 + A^°-^ +1 a n+1 c ) 



1, Xclquj~ 



-70 



It means that we can imitate different boundary spins by one L-operator 

V u n+ ic n+ i, AaoCn+iu; ,u J 

which is of the form ([6]) with special parameters. Thus in this paper we consider the 
monodromy matrix defined by (TTTT) for the quantum chain with n + 1 sites, where L n+1 (A) 
is given by (jlfip . The corresponding transfer-matrix with fixed boundary spins is 

t B (A) = tr T„, +1 (A) = ( 1, Aw-woo ) L X (A) L 2 (A) • ■ ■ L n (X) ( ^ ) . (17) 

The problem of construction of the eigenvectors of the transfer-matrix t B (A) will be 
solved as follows. First, we find the eigenvectors of -Bi(A). Then the eigenvectors of B m (\) 
are obtained recursively as linear combination of the eigenvectors of _B m _i(A). After 
that the eigenvectors for the transfer-matrix (fTTj) are constructed as linear combinations 
of the eigenvectors of B n+ i(X) (the auxiliary model). The multi- variable coefficients of 
this decomposition admit the separation of variables and can be written as products of 
single-variable functions, each satisfying a Baxter difference equation which can be solved 
explicitly. 



3 Eigenvalues and eigenvectors of B m (X) 
3.1 Eigenvalues of B m (X) 

We start from the problem of construction of the eigenvectors of B m (X). According to 
[B m (X), B m (fjL)] = 0, we are looking for the eigenvectors of B m (X) not depending on A 
with the eigenvalue being a polynomial in A. Factorizing this polynomial we get 

m— 1 

B m (X)* Pm = Ar m , ur^° H (A + r mtS u-^) V Pm , (18) 

s=l 

where r m s , s — 0, . . . , m — 1, is a set of constants and we shall use the phases 

Pm = (Pm,0, Pm,m-l) € (Z N ) m (19) 

as labels of the eigenvectors. These formulas are valid for all m = 1, . . . , n + 1, that is 
including the chain with boundary L-operator (fT6|) . As we shall see later, in the latter case 
not all the phases p n +i are possible but only those which satisfy the following restriction 
related to the boundary spin 7 : 

n 

Pn+i ■= 7^ Pn+i,s = 7o + 1 modiV. (20) 

s=0 
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Let us define the "averaged" counterpart 0(X N ) of a quantum cyclic operator 0(A) 
using averaging procedure [12] 

0(X N ) = (0}(X N )= U s ez N 0(uJ s X) (21) 

and apply this procedure to the entries of the quantum L-operator ([6]). Denote the result 
hyC k (X N ) 

r ixN.f l-en%X N -tX N (a%-b%) \ 
k{ } ~\ $-d» b N k d N k j^ - eX N a%c% J ' {Z2) 

where e = (— 1)^, and call it as the "averaged" ^-operator of the BBS model. In partic- 
ular, the averaging of (1TB]) gives 

C n+ i(X N ) = ( n _ c \A f /7 A r 3v ) • (23) 

V c n+l tA u L n+1 / 

Accordingly, the averaged monodromy T m for the m-site chain is 

T m {X N ) = &{\ N ) C 2 (X N ) ■ ■ ■ £M N ) = ( I^n] ) , (24) 

where the entries are polynomials of A . By Proposition 1.5 from [12] (see also [13J), 
these polynomials coincide with averages (A m ), (B m ), (C m ) and (D m ) of the entries of 
(TTTT) . This proposition provides a tool for finding the iV-th powers of the amplitudes r m s : 
applying (|2T|) to ffTHT) we obtain 

m— 1 

B m (X N ) = (-e) m X N rZ j0 H(X N - e<J . (25) 

8=1 

This relation together with (122]) and (124]) allows to find s in terms of the parameters a k , 
bf?, c k , d k and x k , k = 1, . . . , m. The problem of finding the amplitudes r m>s is reduced 
to the problem of solving a (m — l)-th degree algebraic relation. As shown in [1], in the 
case of the homogeneous BBS chain model the problem is reduced to solving a quadratic 
equation only. The described procedure gives the amplitudes r m>s up to some roots of 
unity. In fact we can fix these phases arbitrarily because this leads just to relabeling of 
the eigenvectors. In what follows we suppose that we fixed some solution {r m s } in terms 
of the parameters a k , b k , c k , d k and x k . Again in the case m = n + lwe can not define 
the phase of r n+10 arbitrarily. We fix it by the relation 

n 

fn+1 ■= Y\_ r n+l,s = Ud . (26) 

Then the relations ([20]) and (126]) provide correct coefficient a^uj" 10 at A (the lowest term) 
in the eigenvalue of B n+1 (X). 

3.2 One-site eigenvectors for the auxiliary problem 

In BBS model a very important role is played by the cyclic function w p {~f) [10] which 
depends on a Z^-variable 7 and on a point p = (x,y) restricted to the Fermat curve 
x N + y N = 1. We define w p (j) by the difference equation 

Wp(0) = 1 ; 7 e %n- 



w p (i — I) l-u~ix 
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The Fermat curve restriction guarantees the cyclic property w p {pf + N) = w p (7)- The 
function w p {^f) is a root of unity analog of the g-gamma function. 

It is convenient to change the bases in the spaces V k . Instead of \^)k, 7 £ we will 
use the vectors 

which are eigenvectors of the upper off-diagonal matrix element of the operator L k (X): 

A u k \a k - b k v k ) ^ = Xr k (28) 

The coordinates of the Fermat curve points p k = (x k ,y k ) are defined as follows. Let us 
fix some value of r k to satisfy r k = a k — b k . Then 

x k =r k /a k , Uk = b k /a k . (29) 



In the case k = n + 1 that is for (1161) . the eigenvectors (127]) degenerate as follows. 
We fix x n+ i = dJ 70_7n+1 , Dn+i = 0, r n+ i = ao- It gives, up to an inessential constant 
multiplier, w Pn+1 (-f) = 5 7i7n+1 _ 70 and therefore V&t? = \Pn+i + 7n+i -7o)n+i- Th e vector 
Vw+i'' satisfies ( 1281) with = n + 1. On the states with fixed eigenvalue u n+ i = u; 7,1+1 we 
have to identify p n+ i = 70. 

Note, if r k = (in particular, in the superintegrable case) it leads to x k = 0, y k = 1. 
In this case the cyclic function w p {^f) degenerates and fT27|) does not give a new basis in 
V k . So we consider the generic parameters such that r k 7^ 0. 

The operator v fc acts as cyclic ladder operator: 

^ = 4 k Xi • (3°) 

Using f[2"8"j) for k — 1 and comparing to (1T8T) . we write one-site eigenvector as ^ P10 '■ = 
"0pi!o- With ri >0 = r\ we have 

£i(A)# plj0 = A n, a;"' 1 - * Pli0 , v4i(A)^ Pl , = tf^,, + Axi* pi 0+1 . (31) 

3.3 Fermat curve points in the formulas for the eigenvectors of 

B m (X) 

The formula for the eigenvectors of B m {X) is defined in terms of w p {^) function depending 
on four types of points on the Fermat curve x N +y N — 1: 

Pin (^mjl/ro); Prn,s {%m,si 2/m,s)j Pm,s (5m,s; |/m,s)] Pm' ,s' (p^m'.s' ' Vm' s') ' (^2) 

The coordinates of these points are expressed in the terms of amplitudes r mjS , m = 1, . . . , 
n+ 1, s = 0, . . . , m — 1 (defined as some solutions of equations (125]) . m — 1, . . . , n+ 1) by 

~^m',s' Tm,s / Tm' ,s' i •^m,s ^m^rnX m,s / '^mi *^m,s ^m/(^m''m''m,s)) S, S ^ 1. (33) 

In the construction of the eigenvectors of B m (X) we need the points of type p m>s and p m , s , 
s > 1, only for m < n. For all the other types of points m runs up to n + 1. 
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The values of y™f g ,, y m>s , y m>s are defined by the condition on p™f s ,, p m ,s, Pm,s to 
belong to the Fermat curve and the following relations (1 < I < m — 2) on phases of 
y-coordinates: 

_ m— 2 m—l,l rrtn-1 m,k 

l,s llfc=l Vm-1,1 _ 1 / Q/1 x 
c; „ v h s> n J. J. m— l.s rrm-3 m— li ' V ' 

r m _ 2 r m _i j0 r m m _i c m _i ?/ m „i,i y m _i,i ^ y m _ 1 ' J || s=1 y m _ 2 ; s 

where 

'"m ^m,0 ^m,l • • • Tm,m— 1 ■ (35) 

Practically we may define y™f s ,, y m ,s arbitrarily and the coordinates y TO)S by (154"|) . Then 
the points p m)S , due to (|25|) . will belong to the Fermat curve automatically [4J. 
The coordinates of the points p m fl and p m , 1 < m < n, are defined by 

%m,Q r m,0 = r m-l,O a mCm, 2/m,0 r m,0 = ^1^2 ' ' ' K m~\ r m j (36) 

^mXm ''mi 2/m^m ^rn^mTm— l/^w • (37) 

We use the same formulas for the case m = n + 1. Due to (1261) and r n+ i = ao, the 
solution of (157)1 is x n +i = uj^ 1 , y n +i = 0. Correct limiting procedure gives (wp m (p))~ 1 = 
$p,o- We define the coordinates of the Fermat curve point p n +i,o by 

^n+l,O r n+l,0 = r n,0 a 0C n +l, 2/n+l,0 r n+l,0 = Xl><2 ' ' " Xn«0 • (38) 

Note that from (136|) we had to define in+^o by relation avt-i,o?Vi-i,o = ^n,o a oCn+i^ 7n+1 ~ 70 
but we want to avoid entering the values of boundary spins into the definition of Fermat 
curve points. Such change of definition of x n +i,o means that instead of W( Xn+1 0tVn+1 )(p) 
we have to use 

/ \ _ ^Pn+l.o(P + 7n+l -7o) , QO X 

7 " +1 " 7 °— " ^ + ,o(7. + i-7o) (39) 

with point p n +i,o defined by This formula allows to move the dependence on the 

values of boundary spins 7 and 7 n+ i to the argument of W( Xn+1 0) y n+1 )(p)- 



3.4 Recursive formulas for the eigenvectors of B m {\) 

Recall from (131)) that the vector pi := ipp{ i0 G V\ is eigenvector for -Bi(A) and from (1TB]) . 
(TTOT) that the eigenvectors \E r Pm of B m (X) were labeled by the vector p m = (p m ,o, • • • , p m ,m-i) £ 
(Z N ) m . Let us further define: 

Pm — ZlfcLo 1 P™,&> Pm = SfcLl 1 Pm,fc; P m = (Pm,l, • • • ,Pm,m-l) G (2W) m (40) 

The vector denotes the vector p m in which p m ^ is replaced by p m) k ± 1, i-e. 

Pm = (Pm,0, • • • , Pm,fc ± 1, • • • , Pm,m-l), k = 0, 1 , . . . , 771 - 1 . 

The following Theorem [TJ proved in [1] gives a recursive procedure to obtain the eigen- 
vectors \l/p m G )A m ), 2 < m < n, of B m (X) from the eigenvectors $ Pm _; G l/™ -1 ) of 
-B m -i(A) and single site vectors ip^ G V m defined by (l2Tj) . The recursion is starting 
from the already defined ^ Pl0 - The Theorem 1 is valid provided ^ 0, the polynomi- 
als B m (\ N )/ X N , m = 2, ... ,n, have nonzero simple zeros and det7^(er^J ^ (cf. the 
definition of the 5-representation in [T2]). 
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Theorem 1 The vector 

*P™ = Q{Pm-liPm\Pm)y pm -,®^ (41) 



Pm _iS(z JV ) m - 1 

Pm £^jV 



where 



^j(Pm—f>m-l){pm—pmft) 

Q\Pm-\i Pm\Pm) = 7 7\ Tz rr X 

W Prn0 {Pm,0 - Pm-1,0 - L)Wp m (p m ~ Pm ~ 1) 
rm— 2 t r 772 — 1 



X 



11™=? (Pm-lJ - Pm-u) f_7 w p m -i,i(Pm-l,0 

is eigenvector of B m (X): 

m—l 

B m {\)^ Pm = Xr mfi uJ- p ™' J] ( A + r m,kUJ~ p - k ) # Pm . (43) 

fe=i 

t/ie m — l zeros \ m ,k of the eigenvalue polynomial of B m (X) 

X m ,k = -r mtk uj- p ^\ k = 1, . . . , m - 1, (44) 

the operators A m and D m act as shift operators for the k-th (and 0-th for D m ) index of 

A m (Kn,k) * Pm = MP'J * p +* > D m (X m>k )^ Pm = <p k (p'J ^ pt o,- k , (45) 
where 

MP'J = ~— F m (X m ,kH II C\- > (46) 



s=l 



5' -1 m_1 

r m uj p ™ 



MfL) = n m-2 m,k II F -( A --fc) , (47) 

' m ~ 1 ll s =l Z/m-l.s s =l 



^s(A) = (b s +ua s K s \) (Ac s +4/x s ). (48) 

The next theorem gives the formula for the eigenvectors of B n+ i(X) in the space 
(which is restriction of the initial space of states y( n+1 ) to the subspace of fixed value jn+i 
of spin at (n+l)th site) in terms of the eigenvectors of B n (X). As was explained before, the 
components of p n +i are not independent but satisfy (1201) . So we will label the eigenvectors 
#5 of B n+1 (X) by the set p' n+1 = {p n +i,u Pn+i,n)- Then p n+ i i0 = 70 + 1 - p' n+v 

" n+1 

Theorem 2 The vector 



*JL = E Q B (Pn|PUl)^P. (49) 



-"n + l 



/ Ct )('W)+ 1 -Pn)0 / n+l- 1 ) 

Q (Pn\p' n +l) ^ 



X 



Wpn+i,o(7n+l-^H-i-Pn,o) 
Ul=im=^ p: yAPn,l-Pn + l,k) Wpn j (-P„,) (5Q) 
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is eigenvector of B n+ i(X): 

n 

B n+1 (X)^ B p , n+i = Xr n+lfi u-^° H (A + r n+ljk u-^) # B ;+i . (51) 

k=l 

At the n zeros X n +i,k of the eigenvalue polynomial of B n+ i(X) 

K+i,k = -r n+hk uj~ p "+ 1 - k , k = 1, . . . , n, (52) 

the operators A n+ i and D n+ i act as: 

A n+l (A n+ljfc ) ^ B p , n+i = , D n+1 (A n+ i, fc )*J n+i = ^fc(pUi) • ( 53 ) 

Proof. The formula for the eigenvectors \I/ B , of 5 n+ i(A) follows from the for- 
mulas (I4ip and (142p at m = n + 1 after appropriate limiting procedure: the factor 
w P„+i(Pn+i — Pn+i ~ 1) kills summation over p n+ i (fixing p n+ i = 70) and becomes 1, 
the factor w Pn+10 (p n+lj0 - p n ,o ~ 1) due to ([39]) becomes w Pn+10 (j n+1 - p' n+1 - p nfi ) up to 
an inessential constant multiplier w Pn+1 (7 n +i — 70) • The action formula for A n+ i follows 
from the identity Xa Lj~' yo A n+ i(X) = B n+ i(X). □ 

4 Eigenvalues and eigenvectors of the boundary 
transfer- matrix t B (A) 

After having determined in Theorem 2 the eigenvectors \I/ B / of the auxiliary system, 

we are looking for the eigenvectors of t B (A) from (1TT1) as linear combinations of the 
eigenvectors \I/ B , 

Pn + l 

From the expressions for the //-operators we can see that the eigenvalue of t B (A) is 
t B (A| 7 o, 7n+i, E) = E + E X X + ■ ■ ■ + E n X n + E n+1 X n+1 , (54) 
where E = {Ei, . . . , E n } and the values of E and E n+ \ are 

n 

E = 1, E n+1 = Y[ a m c m ■ a Q c n+ iu ln+1 ~ 10 . (55) 

m=l 

The possible values of E can be found from the functional relations. We define the 
following polynomials depending on unknown set E: r^(A) = 0, r^'(X) = 1, r^ 2 ^(A) = 
t B (X) and recursively 

T ^ +1 \X) = T {2 \u j - l X) r (j) (A) - u p zifJ-^X) r " 1} (A), j = 2,3,...,N, (56) 

where z(X) = JXJii ^.(-V^O- Then the "truncation" identity 

r<" +1 >(A) - u»z(X) r^-^wA) = A n+1 (X N ) + V n+1 (X N ) (57) 

defines possible values of E. Note, the polynomial A n+ i(X N ) + V n+1 (X N ) corresponds to 
a q + a q in [5] . 
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As was mentioned in [5], in the case of parameters considered in the present paper 
we have z(X) = (due to F n+ i(A) = 0) and therefore all the set of functional relations 
reduces to one relation (we use the averaging (I2TI) ): 

(t B )(X N ) = A n+ i(X N ) + T> n+ i(X N ) (58) 

or, equivalently, (A n+ i+D n+ i) = (A n+1 ) + (D n+1 ). At this point it is natural to conjecture 
the formula for the eigenvalues of transfer- matrix [5]: 

n+l , A 



t B (X\<r n+1 ) = H — ^ + 1 , (59) 
i=i J 



where <x n+1 = (ax, cr 2 , . . . , a n+ i) E (Zn) and 

(J\ + a 2 H h a n+1 = 7 n+ i - 7 mod N, (60) 

The numbers Sj, j = 1, 2, . . . , n + 1, satisfy two relations: 



n+l n 

n s J~ i = n am ° m ■ a o c «+i> (6i) 

j=l m=l 

n+l / 

n I 1 - (-v N w) = a ^ x ^ + v ^ n )- (62) 

The relation ( |62l) follows from functional relation ( !58l) and averaged ( l59l) . The relation 
( 1621) means that (— l) N Sj are roots of polynomial A n+ i(X N ) + ©„ + i(A Ar ) with respect to 
A^. Note, the relations ( 1601) and ( 16TI) provide correct value of E n+ \ in ( 1551) . 

In the paper [5], Baxter conjectured that to each set er n+1 satisfying (I60l) there corre- 
sponds the eigenvector ^o-n+i of transfer- matrix t B (A) with eigenvalue (159|) . The method 
of functional relations gives in principle the eigenvalues without information on their mul- 
tiplicities (0,1 or more). The method of separation of variables developed in pE] for BBS 
model gives the formulas for the eigenvectors if the eigenvalues are provided. In particular, 
more precise information on the multiplicities of the eigenvalues can be obtained. 

According to [4], we are looking for Q<y n+1 to be of the form 



= 2^ Q>n + iK+i) *p; +1 , (63) 

P'n+l 

where n 

n B („> l~ ^ rife=l Qk{Pn+l,k) 



rim,m' = l W n+l,m' (p n +l,m ~ Pn+1, 

(mytm') "n+l,m 



and the functions of separated variables q k (p n +i,k) satisfy the Baxter type difference equa- 



tions, k — 1, 2, . . . , n: 



^ B (A„+i,fc|cr n+1 ) q k {p n +i,k) = A^(X ntk ) q k (p n+ltk + l) + A k (ujX n+ltk ) q k (p n+hk -l) (64) 
with 

■n 

AjJ-(A) = Xk 1 (A/w)-" F m (A/w); A, (A) = Xk (X/u) n F n+1 (X/u); (65) 

m=l 
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n n—1 

r n+l,O r n / TT n+l,rre / n+l,k \ TT ra+l,fc lrr\ 

Xk = (11 Vn+llk /yn+llm) [[ Vnjn' ■ (66) 

u ' n+1 



m=i m=l 



Since in the case of parameters considered in this paper we have F n+ i(X) = we get 
A A T(A) = 0, and equations (JM|) can be solved explicitly in terms of functions ^(7). The 
answer for Q B (p' n+1 \cr n+1 ) is 

nB , , , , nLin^Npgg^j-Pn+i,*) " ™ Wp B ( Pw+1|fc -i) 

Q (Pn+lkn+l) = ™ J V II II 7 TT' ^ 6T ) 

llm,m'=l W n+l,m' [Pn+1 ,m ~ Pn+l,m' ) W p B , \~Pn+l,k ~ 1 ) 

(m^vn.1) Pn+l,m fc=l 771=1 ^m,fc 

where the Fermat curve points pj®, and are defined by 

BB / B l // \ ~B // 

•Ejk '*n+l,k/^jj *^m,fc "m/ \P"m%mTn-\-'\.,k ) j •^m,k ' n+X,k^"m^m/ ""m 

and conditions on the discrete phases of y?f , : 

n n 

JJ ua m d m y B j,y B j, = Xk Y[ VjS > fc = 1, 2, . . . , n . 

m=l j'=l 

The validity of these relations up to iVth root of unity is provided by (1621 . Straightforward 
calculation gives 

kB/\^ _ +B 



t B (A)$ <Tri+1 =t B (A|o- n+1 )$ CTll+1 , 



where the eigenvectors $ CTn+1 and corresponding eigenvalues are defined by (j551) . (j57)) . (j59|) . 

5 Ising-like Z^- model with fixed boundary spins 

In this section we fix special values of parameters b k = d k = of L-operators ([6]): 

1 + Ax fc v fc , Au" 1 ^ 
L k (\)=[ ), k = l,2,...,n, (68) 



and consider the transfer-matrix (1171) of the BBS model with fixed boundary spins. The 
coefficients of this transfer-matrix expansion in A give a set of commuting operators. 
In the previous section we gave the formula for the common eigenvectors of this set. 
The coefficient at A gives Hamiltonian of Ising-like Zat quantum chain model with fixed 
boundary spins [6j [TT] 

n n+1 

Hi = x fe v fe + 2j a k c k+ i\i k 1 u k+1 , (69) 

k=l k=0 

where u = cj 70 , u n+1 = uj ln+1 are fixed boundary spins. The corresponding eigenvalue 

n+1 
3=1 
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follows from ( |59|) . the set cr n+1 satisfy ( |60|) . and the amplitudes Sj have to be found 
from (I6T1) and (1621) . Using the explicit expressions for averaged L-operators (|22|) with 
bk = dk = and (|23|) we can rewrite right-hand side of (162]) as 



An +1 (\") + V n+1 {\ N ) = tr T n+1 {\ N ) = tr A(A W ) • • • C n+1 {\ 



where 



4(A 



1 

c 







-A' 



cf +1 



1 -eA 7 ^. 2 



2fc 



1 - eA"/& 



-eA Ar /3| fc+1 




-eA JV /?| fc+1 
—eX (3 2 k+i 



1 o 

with e = (-1)^, /Wi = {a k c k+ i) N/2 , (3 2 k = xf /2 - Finally we get 



trT n+1 (A w ) = (l,- e A^ 2 )(n( 

U=i ^ 



-e\ N [3 2 



2k 







-eA 






where the multipliers in the product over increasing are ordered from the left to the 
right. Rewriting it as determinant of (2n + 2) x (2n + 2) tri-diagonal matrix we get s^, . . ., 
from the eigenvalue problem for (2n + 2) x (2n + 2) bi-diagonal matrix 



B 



( 


&x 


. 








\ 







02 ■ 













& 


. 
















. 





02n+l 




V o 





. 







/ 



which has 2n+2 eigenvalues ±s^ N ^ 2 , j = 1,2,..., n+1. The method of finding eigenvalues 



of Hamiltonian (1691) in terms of eigenvalues of B was given in [TT] . 

Another model with A n (X) as generation function of the Hamiltonians has similar 
boundary conditions and was considered in [TJ]. A wide class of models with integrable 
boundary conditions can be obtained in the reflection equation approach with the use 
of the cyclic L-operators ([6]). So it is interesting to construct the eigenvectors in these 
models. In particular case of the relativistic Toda chain at root of unity, the formulas for 
the eigenvectors were obtained in [15]. 
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